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Work 


A simplified method of matched asymptotic expansions 
has been developed where the common part in composite 
solution is generated as a polynomial in stretched variable 
instead of actually evaluating the same from outer solution 
(see items (i) and (ii) in the enclosed list of 
publications). This methodology has been applied to the 
solution of the exact equations for three dimensional 
atmospheric entry problem. A composite solution is formed in 
terms of an outer solution, an inner solution, and a common 
solution. The outer solution is obtained from 
gravitationally dominant region, whereas the aerodynamical ly 
dominant region contributes to the inner solution. The 
common solution accounts for the overlap between the outer 
and inner regions. In comparison to the previous works, the 
present simplified methodology yields explicit analytical 
expressions for various components of the composite solution 
without resorting to any type of transcendental equations to 
be solved only by numerical methods. (See item (iii) in the 
enclosed list of publications). 

In the next stage, we address the optimal control 
problem arising in the noncoplanar orbital transfer 
employing aeroassist technology. The maneuver involves the 
transfer from high Earth orbit to low Earth orbit with a 
prescribed plane change and at the same time minimization of 
the time integral of the heating rate of the spacecraft. 

With a suitable performance index, we formulate the optimal 
control problem. Using Pontryagin minimum principle, the 
state and costate equations are obtained, leading to the 
nonlinear two-point boundary-value problem. This problem is 
solved numerically by using multiple shooting method (see 
item (iv) in the list of publications). On similar lines, 
the optimal control problem for coplanar orbital transfer is 
also being investigated using multiple shooting method. 
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During the same period, other related research works 
have been carried out and are briefly mentioned below. 

1. An important work in the same period is the final 
preparation of the forthcoming book entitled, "SINGULAR 
PERTURBATION METHODOLOGY IN CONTROL SYSTEMS, authored by Dr. 
D. S. Naidu, and being published under IEE Control 
Engineering Series, by Peter Peregrinus Limited, Stevenage 
Herts, England. This book is scheduled to appear in February 
1988. (See item (v) in the enclosed list of publications). 

2. As an outgrowth of earlier work on singular 
perturbations and time scales in discrete control systems, 
it has been found that to a zeroth order approximation, 
these two approaches yield identical results. (See item 
(vi) in the enclosed list of publications). 

3. Other work is concerned with the preparation of a 
NASA Technical Publication (see item (vii) in the list of 
publications) . 





(i) D. S. Naidu and D. B. Price, "On the method of 
matched asymptotic expansions", SIAM Annual Meeting and 35th 
Anniversary, Denver, CO, October 12-15, 1987. 

*(ii) D. S. Naidu and D. B. Price, "On the method of 
matched asymptotic expansions", accepted for publication in 
Journal of Guidance, Control and Dynamics, 1988. 

*(iii) D. S. Naidu, "There -dimensional atmospheric 
entry problem using method of matched asymptotic 
expansions", Accepted for presentation at 1988 American 
Control Conference, Atlanta, GA, June 14-17, 1988. 

*(iv) D. S. Naidu, "Optimal control of aeroassisted 
noncoplanar orbital transfer vehicles", Draft paper 

(v) D. S. Naidu, "Singular Perturbation Methodology 
in Control Systems", IEE Control Engineering Series, Peter 
Peregrinus Ltd., Stevenage Herts, England, 1988. (in press) 

(vi) D. S. Naidu and D. B. Price, "On singular 
perturbation and time scale approaches in discrete control 
systems", accepted for publication in Journal of Guidance, 
Control and Dynamics, 1988 (in press). 

(vii) D. S. Naidu and D. B. Price, "Singular 
perturbations and time scales in digital flight control 
systems”, NASA Technical Publication, Spacecraft Control 
Branch, Langley Research Center, Hampton (in preparation). 
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ON THE METHOD OF MATCHED ASYMPTOTIC EXPANSIONS 


Dr. D. S. Naidu* 

Associate Research Professor 
Department of Electrical and Computer Engineering 
Old Dominion University 
Norfolk, VA, 23508 

and 

Dr. D. B. Price* 

Assistant Head 

Spacecraft Control Branch, MS 161 
NASA Langley Research Center 
Hampton, VA, 23665-5225 

INTRODUCTION 

Singular perturbation problems, where suppression of a small 
parameter affects the order of the problems, have been solved by a 
wide variety of techniques* -5 Two of these techniques, singular 
perturbation method (SPM) 1 ' 5 and the method of matched asymptotic 
expansions (MAE) 2 ' 9 have been independently developed to a 

reasonable level of satisfaction. Essentially, the SPM consists 
of expressing the total solution in terms of an outer solution, an 
inner solution, and an intermediate solution. On the other hand, 
in the method of MAE, a composite solution is constructed as the 
outer solution, the inner solution and a common solution. 

In this Engineering Note, a critical examination of the method of 
matched asymptotic expansions reveals that the various terms of 
the common solution of MAE can be generated as polynomials in 
stretched variable without actually solving for them from the 
outer solution as is done presently. This also shows that the 
common solution of the method of MAE and the intermediate solution 
of the SPM are the same and hence that these methods give 
identical results for a certain class of problems. An illustrative 
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example is given. 

METHOD OF MATCHED ASYMPTOTIC EXPANSIONS 
The method of matched asymptotic expansions has been 
extensively used in fluid mechanics.* In this method, a 
composite solution is expressed as an outer solution, plus an 
inner solution, and minus a common solution. 


We describe briefly the method of MAE as applicable to 

initial value problems. Consider 

= f ( x , z, £, t) (la) 

*dt = g(x ’ z> £ ’ (lb) 

where x, and z are n~and m-dimensional state vectors respectively 
and c is a small positive parameter responsible for singular 
perturbation. We begin by representing the solutions in the form 
of a series in powers of £ as 

OD CD 

x(t,«) = E x u> (t)« v ; 2 (t,c) = E z‘ l> (t)f l (2) 

isO 

and determine the various terms x <u (t) and z <u (t) by means of 
formal substitution of Eq. (2) in Eq. (1) and comparison of 
coefficients of equal powers of £. Then the following set of 
recursive equations are obtained. For zeroth order approximation, 


dsltlm , f . [x .o. (t)> 0> t] 


dt 


a O *► <o> i . * CO) . , i . . « 

o = g [x (t), z (t), o, t] 


(3a) 

(3b) 
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and for first order approximation, we have 


, <*>... 

dt 


. <*>> v 

dz _(t) 
dt" 


= f* [X ci, (t) f X <0, (t), 2 <0> ( t ) , t] (4a) 

= 6 [x 4> (t), z W> (t), x <0> (t), z‘ 0> (t), t] (4b) 


where the notation f°, and f 1 is used to indicate all the terms 
on the right hand aide. Since the series of Eq. (2) corresponds 
to the solution outside the boundary layer, it is called an outer 
series. 

The solution of Eq. (3) is obtained by using x (t=0) = 
x(0); and in general z <o> (t=0) * 2(0). On the other hand, the 
solution of Eq. (4) poses a problem, since the initial condition 
x w> (t=0) is not yet known. Once x a> (t) is solved for, 2 W> (t) 
is automatically known from Eq. (4b). In order to relate the 
outer series of Eq. (2) to the solution of (1) in the boundary 
layer, we use a stretching transformation 

r = t/c (5) 

Then using Eq. (6) in Eq. (1), the stretched or inner problem 
becomes 

ffLt) = f [x(t ) , z(r), s, *t] (6a) 

= g[x(r), z(r), e, ct] (6b) 

This has inner series expansions of the form 

x(t,£) = £?’(T)r l ; z(t,c) = £z <x> (t)c 1 (7) 

isO 1 = 0 

Substitution of Eq. (7) in Eq. (6) and comparison of coefficients 
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result in for zeroth order approximation as 


,- < o > . . 

dr 


0 


(8a) 


.-<0> , X 

dz__.il) 

dr 


— o r -<o> . -<o>, x _ 

= 8 [x (r). z (T)] 


(8b) 


and similarly we get equations for first order approximation. 

The inner problem of Eq. (6) has initial conditions as 

x o> (t=0) = x(t=0); z o> (t=0) = z(t=0 ) (9a) 

x 1> (t=0) - 0; z°(t= 0) = 0; i > 0 (9b) 

Still, we have not resolved the problem of determining the 
initial value x a> (t=0) of the outer problem of Eq. (4). This is 
done by using a matching principle of the method of MAE? ' 9 Thus the 
matching principle is stated as 

inner expansion of outer solution, (x°) 1 = 

outer expansion of inner solution, (x 1 ) 0 (10) 

To any order approximation, the composite solution x is 

C 

given by 

x c = x° + x - (x°) i = x° + x l - (x l )° (11) 

where x , and x are the outer and inner solutions respectively 
to any order of approximation and (x°) 1 = (x l )° is also called 
the common solution. Similar expressions can be given for z also. 

AN EXAMINATION OF COMMON SOLUTION 

In this section, we will show that the common solution 
defined as the inner expansion of the outer solution is simply 
formulated as a polynomial in the stretched variable. The steps 
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involved in obtaining the common solution are (i) express the 
outer solution in the inner variable t, ( i i ) expand it around e - 
0, and (iii) rearrange the resulting solution in powers of e. 
Thus, consider the outer solution as 


x°(t) = x fo> (t) ♦ *x <4> (t) + (12) 

We express this outer solution in the inner variable r = t/c as 

x°(rr) = x°\ct) + ex t} {cT) + (13) 

Expanding Eq. (14) around e - 0, we get 


/ o 

(X ) = 


[ 


<o> 


A < ° > , X 

«->l + ‘-k - L - v - 1 

* £=a 1 £=( 




'[«*"(«) i ♦ | * . . .] 

L c-o ' e-o J 


(14) 


Now evaluation of function x (*t ) at e - 0 in t- plane is 
the same as its evaluation at t = 0 in t-plane, and the partial 
derivative of function x°(*r), with respect to e in T-plane is 
the same as its partial derivative w.r.t. t multiplied by r in 
t-plane. Thus, 


uv = [ x o> (t=0) + +•••] + 

x <4> (t=0) +rT-"— +...J 

= x <0> ( 0 ) + x‘ 4> (t=0) + rJt <O> (0)j +. 


, ° %*■ 

(x ) = 


<o> . . 

X (t ) 


+ *X 4> (t) 


(15) 
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1- ~<0> / V CO) , ,, . "ci) . , <*>._, ,<o) 

where, x (t ) = x (0); x (t ) = x (0) + r* (0), 

and the dot over x denotes differentiation of x w. r. t. t . 

Similar expression can be obtained for the function z. Let us 

note that the intermediate solution of SPM is obtained by (i) 

expanding the outer solution around t = 0, (ii) expressing it in 

the inner variable t, and (iii) rearranging the resulting 

solution in powers of * V® Then, the common solution of Eq. (16) 

of the method of MAE is found to be the same as the intermediate 

solution of the SPM. Thus, the outer and inner solutions being 

the same in the SPM and the method of MAE, we clearly see that 

these two methods give identical results. Essentially, this 

equivalence means that the expansion of the outer solution around 

t = 0 and transformation into t -plane is the same as 

transformation of the outer solution into r -plane first and then 

expansion around e = 0. The main advantage of the present 

formulation of the common solution is that its various terms can 

be very easily generated as polynomials in t and hence one need 

not have explicit outer solution to arrive at the common solution. 

In this way, we suggest an improved method of MAE, where the 
outer and inner solutions are obtained as before and the common 
solution is generated simply as a polynomial in the stretched 
variable r, instead of evaluating it from the explicit outer 
solutions as is done usually* 

EXAMPLE 

Consider a simple second order system so that we can get 
explicit expressions for the solutions. 
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dx _ 
dt z 

x(t=0) = a 

(16a) 

dz 

r dt = - x ' z 

z(t=0) = b 

(16b) 


Applying the method of MAE described in Section 2, we summarize 
the results as follows. The outer solutions corresponding to Eqs 
( 3 ) and ( 4 ) are 

x <0> (t) = ae z <0> (t) = -ae 1 

x (i> (t) = [x (1> (0) - at]e" 1 • (17) 

z <4> (t) = [-x cl> (0) + at - a]e _t 

The inner solutions corresponding to Eqs. (8) and (9) are 

x <0> (t) = a; z <0> (t) = -a + (a + b)e T 
i“’(T) = (a + b) -aT - (a + b)e T 
z <4> (t) = -(2a + b) + aT + [2a + b + (a + b)T]e 

Considering the two-term expansions only, the common solution (CS) 
for x is obtained as 

(CS) x = (x l )° = (x°) 1 (19) 

From Eq. (18), we obtain (x x )°, the outer expansion of the inner 
solution by first expressing the inner solution in the outer 
variable t = er and then expanding it around c - 0 . Thus 

(x l )° = a( 1 - t) + *(a + b) (20) 

Next, from Eq. (17), we obtain (x°)*', the inner expansion of the 
outer solution as 

(x°) 1 = a( 1 - t) + <?x <1> (0 ) (2D 
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Alternatively, in the present approach, we formulate (x°) 1 as 


(xV = x O> (t=0) + * [x w> (t=0 ) + Tk <o> (t=0)] 

= a( 1 - t) + *x tt, (0) (22) 

Equating Eqs. (20) and (21), we get the value of undetermined 
coefficient x u> (0) = (a + b). Similarly for z, we have 

(CS)^ = (z l )° = (z l )° (23) 

From Eq. (18), we obtain (z v )°, the outer expansion of the inner 
solution as 

(z')° = -a(l - t) + r[-(2a + b)] (24) 

Next, we obtain (z°) v , the inner expansion of the outer 
solution as 

(z°) 1 = -a(l - t) + e z <4> ( 0 ) (25) 

Alternatively, in the improved method, we formulate (z°) 1 as 

(z ) = z (t=0) + ^[z (t=0) + rfc (t=0)] 

= -a( 1 - t) + cz {± \0) (26) 

Using Eqs. (23)- (25), we get the value of the undetermined 
coefficient z <4> (0) as 

z <4> ( 0 ) = -(2a + b) (27) 

The composite solution corresponding to Eq. (12) is 

X c ( t , <■ ) = ae _t + r[(a + b)(e _t - e~ l ' c ) - ate" 1 ] (28a) 

z (t,r) = -ae~ l + (a + b)(l + t)e -1/r ♦ 

C 
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CONCLUSION 


In this paper, a critical examination of the method of matched 
asymptotic expansion have revealed that the terms of the common 
solution could be generated as polynomials in stretched variable 
without actually solving for them as it is done presently. We have 
also seen that the common solution of the method of matched 

asymptotic expansion is the same as the intermediate solution of 
the singular perturbation metho^find hence these two methods give 
identical results. Two examples have been given for 

illustration . 
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Abstract: The analysis of a three-dimensional atmospheric entry 
problem using the method of matched asymptotic expansions is 
considered. A composite solution is formed in terms of an outer 
solution, an inner solution and a common solution. The outer 
solution is obtained from gravitationally dominant region, whereas 
the aerodynamically dominant region contributes to the inner 
solution. The common solution accounts for the overlap between the 
outer and inner regions. In comparison to the previous works, the 
present simplified methodology yields explicit analytical 
expressions for various components of the composite solution 
without resorting to any type of transcendental equations to be 
solved only by numerical methods. 
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Nomencl at ore 


: constants of integration for outer solution 

C. : constants of integration for inner solution 

C : drag coefficient 

D 

C : lift coefficient 

Li 

D : drag force 

g : gravitational acceleration 

g s : gravitational acceleration at surface level 

h : nondimensional altitude 

I : inclination of the plane of the osculating orbit 

L : lift force 

m : vehicle mass 

r distance from vehicle center of gravity to planet center 

r : distance from vehicle center of gravity to surface level 

B 

S : aerodynamic reference area 

t : time 

V : velocity 

v : nondimensional velocity 

a : angle between the line of the ascending node 

and the position vector 
ft : inverse atmospheric scale height 

V flight path angle 

V : heading angle 

o : bank angle 

€> : down range angle or longitude 

<p : cross range angle or latitude 

P : density 

p gravitational constant of Earth 

O : longitude of the ascending node 

1. Introduction 

In space transportation system, the concept of aeroassisted 
orbital transfer opens new mission opportunities, especially with 
regard to the initiation of a permanent space station. The 
atmospheric entry problem is of paramount importance for 
aeroassited orbital transfer vehicles (AOTV). 

The atmospheric entry problem involves, in general, the 
solution of nonlinear differential equations by resorting to 
numerical integration. Analytical solutions of a simplified entry 
problem are important from the point of view of serving as a basis 
for investigating more complicated cases and providing a general 
understanding of the structure of solutions. Analytical solutions 
also provide a better foundation for the solution of guidance 
problems. With this in view, attempts have been made to obtain 
approximate analytical solutions for the entry problem using 
asymptotic methods such as the method of matched asymptotic 
expansions, singular perturbation method, and multiple scale 
method [1-11]. Most of these solutions were obtained either for 
the two-dimensional case, or under restrictive assumptions. For 
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instance, in the three-dimensional atmospheric entry problem [9], 
the use of directly matched asymptotic expansions leads to a set 
of transcendental equations which can only be solved by resorting 
to numerical methods. 

In this paper, we address a three-dimensional atmospheric 
entry problem, to be analyzed by the method of matched asymptotic 
expansions (MAE). The solution is expressed in three parts; 
outer, inner, and common solutions. The outer solution is valid in 
the region where gravity is predominant. On the other hand, the 
aerodynamical ly predominant region gives an inner solution. Since 
these two regions are bound to overlap, a matching process is 
required to identify the common solution. Thus, a composite 
solution, valid in the entire region, is constructed as the sum of 
the outer solution and inner solution from which we need to 
subtract the common solution. The matching principle, in other 
words, ties the constants of integration associated with the outer 
and inner solutions with given auxiliary conditions. Compared to 
the earlier work [9], the present method has the following 
features: (i) Analytical expressions have been obtained explicitly 
for the outer, inner and common solutions without facing a set of 
transcendental equations which can only be solved by numerical 
methods, (ii) The composite solution satisfies the given auxiliary 
conditions asymptotically . (iii) The common solution can be 
generated as a polynomial in the stretched variable without 
actually solving for it from the inner limit of the outer solution 
or the outer limit of the inner solution. 

2. Equations of Motion 

Consider a vehicle with constant point mass m, moving about a 
nonrotating spherical planet. The atmosphere surrounding the 
planet is assumed to be at rest, and the central gravitational 
field obeys the usual inverse square law. The equations of motion 
for three dimensional flight of the lifting vehicle are given by 
(Fig. 1) [7,9,10,12], 


dr 

dt ' 

Vsiny 

(la) 

d e _ 
dt 

Vcosycosv' 

rcos<£ 

(lb) 

d<ft _ 
dt 

Vcosysinv / 

r 

(lc) 

dV 

dt " 

- £ - gsira- 

(Id) 

rd r _ 

dt ~ 

LCOSO' „2 , ^ 

m (g V /rjcos r 

( le) 
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(If) 


„dv ( _ Lsino V 
v 3t - ~ jcosj-cosytan# 

where it is assumed that the aerodynamic drag and lift are 


D = |pSC d V 2 ; 




( 2 ) 


the gravitational field is 


2 

g = — — = P/r ; 


P = g r 


(3) 


and the atmosphere is given by 
P = P s exp[-7?(r-r )] 


(4) 


For any particular flight program, the control functions C , 

La 

C D , and ^ are given functions of time and the solution of ( 1 ) 
requires prescribing six initial conditions. 

It is convenient to eliminate time t in (1). Then, we get 


de 

cosy' 

dr 

rcos<j!>tanr 

dj' 

sinv' 

dr 

rtanr 

dV*_ 

psc^v 2 

dr 

msinr 

dr _ 

pSC l cosc' 

dr 

2msin^ 

dy' 

pSC L sino' 

dr 

2msinrcosr 


cosy'tan# 

rtanr 


(5a) 


(5b) 


( 5c ) 


(5d) 


(5e) 


Solution of the set of five first order nonlinear 
differential equations (5) requires integration by numerical 
methods. The aim of the present paper is to obtain approximate 
analytical solutions to (5) using some simplifications in the 
method of MAE. 
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3. Method of Matched -Asymptotic Expans 1 onsC MAE) 

In applying the method of MAE to the three-dimensional entry 
problem, we consider separately the flight in an outer region near 
the vacuum, where the gravity force dominates, and an inner region 
near the planetary surface where the aerodynamic force is 
predominant. There is bound to be an overlap or common region 
where both outer and inner solutions are approximately of equal 
strength. A matching principle is invoked to obtain the common 
part. An approximate solution called the composite solution valid 
over the entire region, is constructed from the outer, inner and 
common solutions. We see that the various components of the 
composite solution are separated depending on the altitude, and 
hence it is appropriate to choose the altitude as an independent 
variable for obtaining the solution of (5). 

Let us define the following dimensionless quantities, 

h = ( r-r )/r ; v = V*/g r ; c - 1/7? r 
6 8 6 6 8 

B = P s SC D /2rV?; X = C L /C D 

Here the constant ftr* is large, i.e., for Earth’s atmosphere ftr* - 
900, and hence the parameter c is a small quantity. 

F rom ( 3 ) - ( 6 ) , we get 

de _ cosycoty 7 

dh " (l+h)cos<£ ' ’ 



d <P 

sinvcoty 


dh ‘ 

( l+h) 


dv 

2Bvexp( -h/c ) 

2 

dh " 

£ siny 

( l+h ) 2 

dr _ 

BXcos<yexp( -h/c ) 

+ F 1 1 Icotr 

dh 

«siny 

( <1+h) va*h/J 

dy 

BXsincyexpt -h/c ) 

cosv'tan^coty 

dh " 

csinrcosr 

(l+h) 


(7b) 

(7c) 

( 7d ) 


(7e) 


Although (7) is ready for analysis by the method of MAE, it is 
more convenient to replace the set of variables (©,tf>,y') with a new 
set of variables (a,n,I) which are related as (Fig. 2) [9,10], 
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COS I = COS0CO8V' 
sin <p - sinlsinot 
COSOt = cos<Acos(9~0) 


(8) 


) 


where I is the inclination of the plane of the osculating orbit, n 
is the longitude of the ascending node, and a is the angle between 
the line of the ascending node and the position vector. 


Using (8) in (7), we get 


da _ _ sinarB^sintyext ^ ~h/£ ) 1 + cotr 
dh “ tanl \ ^sinrcosr J (1+h) 

do _ sine< fB>^sin<yexp( -h/g H 
dh “ sinl \ rsinrcosr J 

dl _ ~ h / g > 1 

dh l rsinrcos r J 


dv 

dh 


2Bvexp( -h/r ) 
rsinr 


2 

(1+h) 2 


dr 

dh 


BXcoso , exp( -h/g ) 
esinr 



1 

(1+h) 


v( 1+h) 2 


cotr 


with initial conditions, a , O , I , v , and r 


(9a) 

(9b) 

(9c) 

(9d) 


( 9e ) 


3.1 Outer CKeplerian) Region 


The outer expansions describe the solution in the region near 
vacuum where the gravitational force is predominant. These are 

assumed as 


a 

0 

1 
v 
Y 


a ( h ) + ccx ( h ) + 

O 1 

O o (h) + rO t (h) + 
I (h) + rl (h) + 

O 1 

v ( h ) + cv ( h ) + 

O 1 

r (h) + cy (h) + 

O 1 




( 10 ) 


By substituting the outer expansions (10) into the original set of 

equations (9), and equating coefficients of c° on either side, the 
set of equations for zeroth-order approximation is 


6 



do© 

coty^ 


dh~~ = 

(1+h; 

(lla) 

dn 



© 

he - ~ 

0 

(lib) 

dl 



© 

dh~ = 

0 

(11c) 

dv 

O 


© _ 
dh 

Lt 

( 1+h) 2 

(lid) 

d ^c 

r i i l 


dh~ = 

L (1+h) v 0 (1+h)Z J co r ° 

(He) 


Let us note that the zeroth-order equations (11) are alternatively 
obtained by letting the small parameter e tend to zero in (9). The 
effect of making £ (=1 /ftr ) zero is that the atmospheric density p 

{= P s ex p(-h/r)} becomes zero and the resulting equations (11) 

describe the region near the vacuum. 

Solving (11), 


v o " 2 [ c i + (l+h)] 
c 2 

COS Y = 

° ( l+h)4v"~ 

o 


C*/( l+h)-l 

cos (a -C_) = 

O 9 ' 

J 1+2C C* 
1 2 


n 


(12a) 

(12b) 

(12c) 
( 12d) 


I 



( 12e ) 


where Q are the constants of integration to be determined. The 

first and higher order solutions are all equal to zero because at 
high altitude, in the limit, the atmospheric density is zero and 
the motion is Keplerian. 
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3.2 Inner C Aerodynamic!) Region 


The inner expansions are introduced to study the limiting 
condition of the solution near the planetary surface where the 
aerodynamic force is predominant. These are obtained by first 
using a stretching transformation 

h = h/c (13) 

in (9) and then taking the limit c •* 0 . This corresponds to the 
region near h = 0, i.e., planetary surface. Thus the stretched 

system becomes 


da _ 

sina fB-\sxno'exp( -h)l 

, ccoty 

dh 

tanl\ sinrcosy J 

( 1 +£*h ) 

do 

sine* fBXsin<yexp( -h)"| 


dh 

sinit. sinycos?- / 


d| 

co - 5 P'Si n^exp ( - h ) ' 1 


dh 

LOSw 1 _ 1 

\ sinreosr } 



(14a) 


(14b) 


(14c) 


dv 2Bvexp( -h ) _ 2c 

dh sin^ (1+ch) 2 


( 14d) 


d r _ B\cos?exp_(-h) + f 1 1 | CJt ~ 

dh sin?' \(l+ch) v(l+ch)J 


( 14e) 


Let us note that x and o% assumed to be external control inputs 
do not undergo transformation. Let the inner solution be expressed 
as 


a 

0 

1 

v 

r 


a ( h ) + ra ( h ) + 

© i 

n (h) + c0 4 (h) + 

I (h) + ci (h) + 

© 1 

v (h) + tv (h) + 
© 1 

r (h) + cy (h) + 
© 1 


(15) 


As before, substitution of (15) into (13) and collection of 

coefficients of powers of c° on either side gives the zeroth-order 
approximation as 
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da 

dh 

dO 

( 

dh 

df 

< 

dh 


sina fE\sinaexp( -h) 


ia fE 1 
© 

* 


tanl i sin r cosy 

© o © 


) 


sin« fB^ 
sinl t s 


sinaexp( -h) 


in r cos r 

O' o © 


) 


fB^sinaexp( -h) 

cosa 


f 


sinr cos?' 

O O 


) 


dv 2Bv exp(-h) 

O O 

dh sinr 

O 

dr BXcosaexp( -h) 

O __ 

dh sin^ 

O 


Solving (16), we get 


(16a) 


( 16 b) 


(16c) 


( 16d) 


( 16e) 


V o = C ± exp (~2r c /^ cosa) 
cosr o = B^cosaexp( -h) + C 2 
sina sinl = sinC 

<5 0 9 

cosa = cosC cosfC -O ) 

O S '4 0'' 


cosl^ = cosC 9 cos|tanalog[tan(r c> /2+rT/4 )] +C 5 | 


(17a) 

(17b) 

(17c) 

(17d) 

(17e) 


where C^are the constants of integration. Here, the sequence of 
solutions is f from (17b), v from (17a), I from (17e), a from 

o o o o 

(17c), and finally O from (17d). We now have the zeroth-order 

O 

outer solutions (12) with C i as the constants of integration and 

the zeroth-order inner solutions (17), where C are the 

integration constants. These constants are determined by a 
matching principle. 

4. Matching Principle! 1 3-1 5J 

Matching is based on the notion that the outer solution valid 
in the Keplerian region and the inner solution valid near the 
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planet surface, must both be valid in some overlap region^. Thus 
matching is accomplished by extending the outer solution into 
the inner region by_ transforming the outer variable h to that of 
the inner variable h (=hA) and taking the limit as c •* 0. This 
is called the inner limit of the outer solution or expansion. 
Similarly, the outer limit of the inner solution or expansion is 
obtained by extending the inner solution into the outer region by 
transforming the inner variable h to that of the outer variable h 
(= eh) and taking the limit as e ■* 0 . By equating the inner limit 
of outer expansion with the outer limit of inner expansion, we can 
determine the constants of integration and hence the common 
solution. A composite solution is formed as the sum of outer and 
inner solutions from which the common solution is subtracted. 

In the earlier work [9], the matching principle yielded a 
relation for the constants <1 in terms of the constants C . Then 

the composite solution is expected to satisfy the given initial 
conditions. This procedure led to the formulation of a set of 
transcendental equations which can only be solved by resorting to 
numerical methods. 

In the present method, we simplify the procedure by asking 
the outer solution to satisfy the given initial conditions and the 
matching principle gives the relation between the constants of 
integration [4], Still, the composite solution satisfies the 
given initial conditions asymptotically. We note that in the 
simplified procedure, we are not faced with any kind of 
transcendental equations and explicit solutions are obtained for 
the composite solution. Moreover, we also get the common solution 
very easily by formulating or generating the various terms of the 
inner limit of the outer expansion as a polynomial in the 
stretched variable as [14,15], 


v m = v (h=0 ) 

© O 

rZ = rjh= 0) 

o o 

« m = C« (h=0 ) 
© o 

n m = n (h=0 ) 
© © 

I* = I (h=0 ) 

o o 


+ £[v 4 (h=0) 
+ r[r i (h=0) 
+ «[« 4 (h=0) 
+ c [O^ (h = 0 ) 
+ «[I 1 (h=0) 


+ hv (h=0 ) ] + 

O 

+ h? (h=0 ) ] + 

o 

+ h*s (h=0 ) ] + 
© 

+ hO (h=0) ] + 

O 

+ hi (h=0 ) ] + 
© 


( 18 ) 


j 


Here the dot denotes differentiation with respect to the 
independent variable h. We note that this is also called the 
intermediate solution in singular perturbation methods [14,15]. 


We now force the outer solution (12) to satisfy the given 
initial conditions, v L , y. <*, a, and I. corresponding to h = h. 

This gives us 


C 4 = v./2 - l/( 1+h ) 


(19a) 
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(19b) 


C 2 = cosr ( 1+h )/v. 


C = cx - cos 

9 t 


C 4 = ft ( 19d) 


C 5 = I ( 19e ) 


Thus, we have the relation between the constants of outer solution 
explicitly in terms of the given initial conditions. In applying 
the matching principle, we first find the zeroth-order inner limit 
of the outer expansion (12), as 


cos y L ( 1+lr )v - 1 


Jl+ ( l+h/)-2j ^cos 2 ?^ ( 1+h. J 


(19c) 


( 20 ) 



The outer limit of inner expansion (17) is 


v; = C 4 exp coso') 

(21a) 

cos?'"' = C. 

o 2 

(21b) 

sina r "sinl m = sinC„ 

o o 9 

(21c) 

cosa m = cosC cos ( C -O m ) 

© 9 4 o 

(21d) 

cos!" = cosC 9 cos^tano'log[tan(f^ , /2+Tr/4 )] +C 5 | 

(21e) 


v" = 2(C +1) 

O 1 

cos?'™ = C 2 /<|2(C j +l ) 
cos(a^-C s ) = (C 2 -l)/ Jl+2C C* 
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Matching (21) with (20), we get constants C in terns of the 
constants C as 

t 


c, = 2( c ,-l >w, Xco()0 , 


{dh c ° r ‘[ c ^ 2<c , +1 >j} 


C 2 = C 2 /i2(C i+ l) 

sinC g = sin£c g +cos -{(c i-i)A sinC 5 

C 4 = C 4 + cos ^cosjCg+cos -f 1+2C, c * })/ coaC 1 


(22a) 

(22b) 

(22c) 


(22d) 


C = cos 

5 


1 jcosC 5 /cosC g j - tano , log^tan^rr/4+cos~ 1 [C 2 /J8(C 4 +1 )jjj 


( 22e ) 


We note that C g in (22d) and (22e) are in turn related with the 
constants C via (22c). 

5. Composite Solution 

The composite solution or expansion is obtained as the sum of 
outer solution (12) and inner solution (17) from which the common 
solution (20) or (21) is subtracted. Thus 


a - ct + a 

CO© 


0=0 + O 

coo 


I = I + I 

coo 


m 

V 

or 

£ o 
i > 

) 


m 

^o< 

or 

— TD 
^o 

) 


m 

°o< 

or 

— m 

a 

o 

) 

► 

o m ( 

o 

or 

a m 

o 

) 


I Tn ( 

o 

or 

I m 

o 

) 



(23) 
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In— 'terms of explicit expressions, we have 


v e = 2/(l+h) + C 4 exp(-2p e />>coscr) - 2.0 (24a) 



Q c = C 4 - cos' 1 [cos«/cosC 3 ] (24c) 


°osI c = cosC 9 cos [tano'log(tan(P o /2+n/4 )) +CJ (24d) 



(24e) 

The above composite solution is expressed in terms of h, 
constants , and C i and the states from the inner solution.. The 
C. s are obtained from (22) and the (1 s are obtained directly from 
the initial conditions via (19). The states of inner solution are 
obtained as explicit functions of the C and h via (17). 

We need to check whether the composite solution (24) 
asymptotically satisfies the given initial conditions. Consider 
(24a) along with (17), (19), and (22). We have 
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As e -* 0 , 


V h =V = rrlirT * 2C * = v i <26b ’ 

Similarly, we can show that r , I . « , and O satisfy their 

c c c c 

corresponding initial conditions asymptotically . 

6. Conclusions 

In this paper, we have addressed the solution of a 
three-dimensional atmospheric entry problem via a simplified 
method of matched asymptotic expansions. The solution has been 
expressed in three parts. An outer solution has been obtained 
in the gravitationally dominant region and an inner solution has 
been formed for the aerodynamically stronger region. A common 
solution has been formed as the outer (or inner) limit of the 
inner (or outer) solution. Finally, a composite solution has been 
constructed as the sum of the outer and inner solutions from which 
the common solution has been subtracted. 

The special features of the present method are (i) The 
composite solution has been obtained in a simplified manner in the 
sense that analytical expressions have been obtained explicitly 
for the various components, outer, inner and commmon solutions, 
without resorting to any kind of transcendental equations which 
can only be solved by numerical methods, (ii) At the same time, 
the composite solution has satisfied the given initial conditions 
asymptotically . (iii) The common solution has been obtained in an 
easier manner by formulating or generating the inner limit of the 
outer solution in terms of a polynomial in the stretched variable. 
The numerical experimentation for the present method is under 
progress . 
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OPTIMAL CONTROL OF AEROASSI STEtf NONCOPLANAR 
ORBITAL TRANSFER VEHICLES 
(Rough draft only) 
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Abstract: The optimal control problem arising in noncoplanar 
orbital transfer employing aeroassist technology is addressed. The 
maneuver involves the transfer from high Earth orbit to low Earth 
orbit with a presribed plane change. The performance index is 
chosen to minimize the time integral of the heating rate of the 
spacecraft. Using Pontryagin minimum principle, the state and 
costate differential equations are derived, leading to a nonlinear 
two-point boundary value problem. This problem is solved by using 
multiple shooting method. 



Nomenclature 


C : drag coefficient 

D 

C : zero lift drag coefficient 

DO 

: lift coefficient 

D : drag force 

g : gravitational acceleration 

g s ; gravitational acceleration at surface level 

K : induced drag factor 

L lift force 

m : vehicle mass 

r : distance from vehicle center of gravity to planet center 

r : distance from vehicle center of gravity to surface level 

6 

S aerodynamic reference area 

t : time 

V : velocity 

ft inverse atmospheric scale height 

Y flight path angle 

V* : heading angle 

<y : bank angle 

9 down range angle or longitude 

<p '• cross range angle or latitude 
P : density 

M i gravitational constant of Earth 
1. Introduction 

In space transportation system, the concept of aeroassisted 
orbital transfer opens new mission opportunities, especially with 
regard to the initiation of a permanent space station [1]. 
In a synergetic maneuver for aeroassited orbital transfer vehicles 
(AOTV’s), the basic idea is to employ a hybrid combination of 
propulsive maneuvers in space and aerodynamic maneuvers in 
sensible atmosphere. Within the atmosphere, the trajectory control 
is achieved by means of lift and bank angle modulations [2-7]. In 
a typical maneuver, we start with a tangential propulsive burn, 
having a characterstic velocity AV d for deorbitting from the high 

Earth orbit and entering into an elliptical transfer orbit. At 

point E the spacecraft enters the sensible atmosphere. The vehicle 

undergoes reduction in velocity due to atmospheric drag and in 

addition, the necessary plane change is performed. At point B, the 

spacecraft leaves the atmosphere augmented by a propulsive burn 

imparting AV for boosting. Once again, the transfer orbit is 
© 

elliptical with a corresponding apogee. Finally, the maneuver ends 


2 


with a reorbit burn having characteristic velocity AV. to make the 

b 

vehicle enter into the low Earth orbit. Thus, the maneuver 
consists of three impulses AV. for deorbit, AV for boost, and av 

do r 

for reorbit. 

In this paper, we address the optimal control problem arising 
in noncoplanar orbital transfer employing aeroassist technology. 
The maneuver involves the transfer from high Earth orbit (HEO) to 
low Earth orbit (LEO) with a presribed plane change and at the 
same time minimization of the time integral of the heating rate of 
the spacecraft. With a suitable performance index, we 
f ormulateoptimal control problem. Using Pontryagin minimum 

principle, the state and costate differential equations are 
derived, leading to a nonlinear two-point boundary value problem 
(TPBVP). This problem is solved by using multiple shooting method 
[ 8 - 10 ]. 

2. Equations of Motion 

Consider a vehicle with constant point mass m, moving about a 
nonrotating spherical planet. The atmosphere surrounding the 
planet is assumed to be at rest, and the central gravitational 
field obeys the usual inverse square law. The equations of motion 
for three dimensional flight of the lifting vehicle are given by 
(Fig. 1), 

$£ = Vsiny (la) 


SI = - ACpV 2 exp( -bft) - (fj/ r 2 )sinr 

= AC L Vcoso , exp( -h/5) + [V/r - /u/r z V]cos r 
= (V/r) cosy's inv' 

= AC L Vsino , exp( -h/?)/cosy - ( V/r )cosycosv'tan<£ 


(lb) 

(lc) 

(l d) 

(le) 


where A = O.SSp /m , h = r-r , p - p exp(-ltf3) and 
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C D = C Do + KC Z for a drag polar. 

3. Performance Index and Hamiltonian. 

For optimal control problem regarding heating, it is required 
to choose the performance index to minimize the time integral of 
the heating rate at a particular point of the spacecraft, say, the 
stagnation point. Thus, the performance index is given by 
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( 2 ) 


J = KjV£ V*‘ 


oa 


dt 


where 1^ is a proportional constant. 

The first step in the optimization procedure using Pontryagin 
principle is to formulate Hamiltonian as 


H = V 3 ' 08 + x^Vsin r 


+ X ■{- AC D V 2 exp( -bft) ~ (M/r 2 )sin/| 


+ X^AC^Vcoso'exp ( ~W ) + [V/r - M/r 2 V]cosyJ 
+ X^{(V/r)cosysiny} 

+ X^AC^ Vsincexp ( -IV? ) /cosy - (V/r ) cosy cosvtan^J 


(3) 


where X's are the costates corresponding to the five states, x’s, 
4. Optimal Controls 

The optimal control equations are given by 


£H . 0 . £H s Q 

ac u ’ 9<y v 

L 

for lift and bank angle leading to 


(4) 


C L = - o>/2KX v V; 


where 


tano' = x^^ycos r 


(5) 


“ = + ( x v / cos? ') 


( 6 ) 
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5. Costata Equations 

The costate (auxiliary) variables X’s are given by 


dX _ _ *H 
dt *x 


( 7 ) 


leading to the corresponding five differential equations as 


dx 


(/?/2)K h V^ # V*°*exp{-^(r-r a )/2} 


- \|(7AC D V 2 exp{-(?(r-r # )} + (2p/r a )sinyj 


* 


+ x^-j(3AC L Vcosyexp{-(5(r-r. )} + (V/r 2 - 2p/r a V)cosy 


} 


+ X^{( V/r 2 ) cosy sinv] + X^|/3AC L V(sino , /cosy )exp{-(5(r-r s ) } 


- ( V/r 2 ) cosy cos^tan^ 


(8a) 


* \ 

dt 


= - 3.08EC h V 2 - oe Vp^ exp{-(J(r-r a )/2} 


X r siny + X y 2AC D Vexp{-/5(r-r»)} 


.{ AC .< 


2 2 

X^-{AC. coso'exp{-(3(r-r_)} + [1/r + #J/r V ]cosy 


} 


>{< 


X -I ( l/r Jcosysinv' 

<P 


} 


- \K 


(sino'/cosr )exp{-/*(r-r )} - ( 1/r )cos^cosy'tan<£f (8b) 


} 


dx 

dt 


y _ 


- x^vcosy + )cosy 
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+ Xy|(V/r - \a/ r 2 V)sinrJ + x^( V/r )ainrsinw 
- xLcj ( tany/cos?" ) sino , exp{ -ft(r~r m ) } 



dX 

dT" 


+ ( V/r )sinT'cos¥'tan£| 

X v/ (V/r)cos^'cosv'/(cos *) 2 

- X ( V/r ) cosycosv' - X (V/r) cosrs inv'tan^ 
<P w 


(8c) 

(8d) 

(8d) 


6. Numerical Data 


The various numerical values used for simulation purposes are 
given below [3]. 

C =0.01; K = 1.11; m/S = 275.0 kg/m 2 

DO 

P e = 1.225 kg/m 3 ; p = 3.986xl0 14 m*/sec 2 

ft = 1/6900 m" 4 ; r £ = 6378.0 KM; = 10.0 

A complete plane change maneuver hs the specifications of 
the initial and final orbit as boundary conditions. The initial 
and final boundary conditions for solving the state and costate 
equations (1) and (8) are given by 


Inital : 

h(t=0) 
V ( t=0 ) 
7(t=0) 
*(t=0) 
V'(t=0) 


= 80 KM 
=7.95 KM/sec 
= -1.25 deg 
= 0 deg 
= 0 deg 


Final : 

h(t=T) 

V(t=T) 

r(t=T) 

*(t=T) 

^(t=T) 


= 80 KM 
=5.0 KM/sec 
= 2 deg 
= 3.83 deg 
= 15 deg 


(9) 


7. Multiple Shooting Method 

The determination of optimal controls (5) requires the 
solution of tenth order, nonlinear, two-point boundary value 
problem (TPBVP) consisting of state equations (1) and costate 
equations (8) and the associated boundary conditions (9). This can 
only be done by numerical methods . The multiple shooting method 
is one of the powerful methods for solving nonlinear TPBVP’ s. The 
corresponding OPTSOL code was developed by DFVLR establishment at 
Oberpf af fenhofen, West Germany. 


In solving any boudary value problem with the given initial 
and final conditions, we asume additional initial data and 
integrate forward so that the solution satisfies the given final 
condition as well. This is also called a simple shooting method. 
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Here, the convergence of the solution is highly sensitive to the 
assumed intial data. It is found that the error due to inaccurate 
intial data can be made orbitrarily small by performing the 
integration over sufficiently smaller subdivided panels within the 
given interval and thereby leading to multiple shooting method. 
Thus, the multiple shooting method is a simultaneous application 
of simple shooting method at several points within the interval of 
integration. Here, the trajectory may be restarted at 
intermediate points using new guesses. Jacobian matrices are 
formed for each segment. The resulting iteration scheme, based on 
reducing all discontinuties at internal grid points to zero, leads 
to a system of linear algebraic equations. 
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